Two-Dimensional Copolymers and Exact Conformal Multifractality 
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Abstract. We consider in two dimensions (2D) the most 
general star-shaped copolymer, mixing random (RW) or self- 
avoiding walks (SAW) with specific mutual avoidance inter- 
actions thereof. Its exact conformal scaling dimensions in the 
plane are derived from an algebraic structure existing on a 
random lattice (2D quantum gravity). The multifractal di- 
mensions r(n) of the harmonic measure of a 2D RW or SAW 
are conformal dimensions of certain star copolymers. The ex- 
act associated f(a) are identical for a RW or a SAW in 2D. 
These are examples of conformal multifractality. 

PACS numbers: 05.20.-y, 02.50.-r, 05.70.Jk, 64.60.Ak 



The concepts of generalized dimensions and associ- 
ated multifractal (Mf) measures have been developed 
more than a decade ago They are encountered 

in many physical situations: strange attractors in dy- 
namical systems, growth phenomena, harmonic measure 
of diffusion-limited aggregates, electron localization, ran- 
dom resistors, and random spin systems. Recently, ana- 
lytic progress has been made for turbulence of a passive 
scalar |i| or for a growth model of DLA || , both where 
the Mf dimensions can be calculated perturbatively. 

It is well known that universal geometrical fractals, 
e.g., random walks, polymers, Ising or percolation mod- 
els are essentially related to standard critical phenomena 
and field theory, for which conformal invariance in two di- 
mensions (2D) has brought a wealth of exact results. By 
contrast, few connections between multifractals and field 
theory have been found, although the algebras of their 
respective correlation functions reveal intriguing similar- 
ities . It remains a challenge to see if an exact descrip- 
tion of some multifractal phenomena could emerge in 2D 
from the conformal invariance classification. 

A particularly interesting multifractal phenomenon 
was singled out some time ago by Cates and Witten 
||. They showed that the moments of the harmonic 
measure, i.e., the Laplacian diffusion field near an ab- 
sorber, the latter taken as a simple random walk (RW, 
i.e., Brownian motion), or self-avoiding walk (SAW, i.e., 
polymer), exhibit in d dimensions multifractal scaling 
for d < 4. The associated exponents can be recast as 
those of star copolymers made of a bunch of independent 
RW's diffusing away from a generic point of the absorber. 
This allowed the perturbative calculation of the Mf spec- 
trum through standard RG theory for polymers and the 
e = 4 — d expansion. 



Star polymers or networks made of self-avoiding walks 
only and their scaling properties are well understood |j) , 
especially in 2D, where all topology dependent exponents 
are known exactly from conformal invariance [p|px|] . For 
simple random walks, scaling exponents describing RW 
non-intersections |Ti]J[2[ ]) conjectured in p3| , have been 
recently derived using conformal invariance and the so- 
called quantum gravity method |p"4j| . 

The aim of this Letter is to extend conformal invari- 
ance to arbitrary mixed copolymers, thereby providing 
the exact Mf harmonic spectrum of a RW or a SAW, an 
example of conformal multifractality. 

Consider a general star copolymer S in the plane 
R 2 (or in Z 2 ), made of an arbitrary mixture of Brow- 
nian paths or RW's (set B) , and polymers or SAW's 
(set V) , all starting at neighboring points. Any pair 
(A,B) of such paths, A, B e B or V, can be con- 
strained in a specific way: either they avoid each other 
(A fl B = 0, noted A A B) , or they are transparent and 
can cross each other (noted AV B), corresponding to four 
different fixed points Jlj|. This notation allows for any 
nested interaction structure; one can decide for instance 
that the branches {Pi € V} e=1 L of an L-star polymer, 
all mutually avoiding, further avoid a bunch of Brownian 
paths {Bk £ B} k=1 n , all transparent to each other: 



(1) 



In 2D the order of the branches of the star copolymer 
does matter and is intrinsic to our (A, V) notation. 

To each specific star copolymer center S is attached 
a conformal scaling operator with a scaling dimension 
x {S) . To obtain proper scaling we consider the Brownian 
paths and the polymers to have the same mean size R. It 
is convenient to define for each star S a grand canonical 
partition function [p| Jll|Jl^] , with fugacities z and z' for 
the total lengths \B\ and \P\ of Brownian or polymer-like 
paths: 



Z R {S) 



E 

B,ves 



z \B\ z r\V\ 



II* (5) 
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where the set of walks of S is constrained by the 
indicatrix lip (S) to stay within a disc of radius R cen- 
tered on the star. At the critical values z c — fig ,z' c = 
fip, where for the RW's /ip is the coordination number 
of the underlying lattice, and yip the effective one for the 
SAW's, Z R decays as ijiT 



1 



Z R (S) ~ R 



-x(S)-x v 



(3) 



where x (S) is associated only with the singularity occur- 
ring at the center of the star where all critical paths meet, 
while x v is the contribution of the independent dangling 
ends or split star. It reads x v — \\B\\ Xb,i + \[P\\ xp t ± — 2V, 
where \\B\\ and \\V\\ are respectively the total numbers 
of Brownian or polymer paths of the star; ig.i or xpi 
are the scaling dimensions of the extremities of a single 
RW x B ,i = 0, or SAW x Pfl = & §,0. The last term 
V = ||Z?|j + |j7 : '|| corresponds in (||) to the integration over 
extremity positions in the disc of radius R. 

When the star is constrained to stay in a half-plane 
with its core placed near the boundary, its partition func- 
tion scales as PJicf 



Z R (S) ~ R 



(4) 



where x (S) is the boundary scaling dimension, x v stay- 
ing the same for star extremities in the bulk. 

Any scaling dimension x in the bulk is twice the confor- 
mal dimension (c.d.) of the corresponding operator, 
while near a boundary (b.c.d.) they are identical: 



2A<°\ 



A(°). 



(5) 



This Letter provides the main lines of a derivation of 
these exponents. 

The idea is to use another representation where the 
RW's or SAW's are on a 2D random lattice, i.e., in the 
presence of 2D quantum gravity |18|. One can indeed 



put any 2D statistical system on a random planar graph, 
thereby obtaining a new critical behavior, corresponding 
to the confluence of the critical point of the infinite ran- 
dom graph with that of the original model. The partition 
function of the random graph G made of, e.g., trivalent 
vertices, reads 



0) = £ 



-0\G\ 



G(x) 



(6) 



where the sum extends over graphs G with a given topol- 
ogy of Euler characteristic x> modulo the group of auto- 
morphisms of G, \G\ being the number of vertices. Near 
the critical point /3 C where \G\ becomes infinite 



Z X (P)~(0-0 D ) 



2-7str(x) . 



(7) 



Tstr (x) = 2 — |x is the string susceptibility exponent 
|l9| . The partition function of the copolymer star S on 
the random lattice G, with the sphere topology (\ = 2) , 
is defined as: 



Z(S)=Y e~^Z G {S) 



(8) 



where the partition function Zq (S) is defined as in (|^), 
with the indicatrix 11r for the star being confined to the 
disc of radius R now being replaced by the indicatrix 1±g 



of the star being embedded in G. The partition sum (g) 
is now three-fold grand canonical, depending implicitly 
on fugacities e~@,z, and z'. A further fugacity term has 
to be added when dealing with a random graph G with 
the disc topology (\ = 1) and boundary length \dG\ , in 
order to define a boundary partition function 



(9) 



where the core star is now on the boundary dG. 



There exists a finite size scaling regime 1 20 where both 
the lattice and the walks become infinite, /?, /?', z, and z' 
approaching together their respective critical values in a 
well defined way. In this regime, the partition functions 
Z, Z, after normalization by the random surface ones 
(0), are expected to scale as pj|: 



\G\ 
\dG\ 



-A(S)-A V 



-A(5) 



\G\ 



-A v 



(10) 

(11) 



is the average size of the ran- 

1/2 

\G\ is the mean length of the 



Z(S)/Z x=2 (0) ■ 
Z(S)/Z X=1 (P) ■ 

Here \G\ ~ (/3 - C )' 
dom lattice, while \dG\ 
boundary. A (S) and A (S) are respectively the bulk and 
boundary conformal dimensions of the star core, dressed 
by gravity. Finally the dimension A v in ( p~0|JTl]) is asso- 
ciated with the star extremities, as was x v in (|3|). Equa- 
tions (0) and ([ll]) are formally identical to and (^) 
in the plane, after recalling (|B|) and identifying the star 
area R 2 in the plane to the random area \G\. 

A general constitutive relation due to Knizhnik et al. 
exists between the conformal dimension A' ' of a scaling 
operator in the plane and the c.d. A of the same operator 
on the random surface A*- ) = A [1 — (1 — A) / n] , where 
k is a parameter related to the central charge of the orig- 
inal statistical model in the plane c = 1 — 6(1 — n) 2 j k 
|p"8| , p2"f . For walks c vanishes, whence k = 3/2, and 

= U (A) , A<°> = U (A) , U (x) = | (1 + 2a?) , (12) 

this relation holding also for boundary operators. 

Here I give a set of basic underlying topological 
"surgery" rules which allow the mixing of geometrical 
operators on a random surface. The set of relations for 
walks is so stringent that it yields immediately the confor- 
mal dimensions of any copolymer star, both on a random 
surface and in the plane. 

Star algebra: The bulk and boundary conformal di- 
mensions, in the presence of gravity, satisfy: 



2A - 7str ( X = 2) = A. 



(13) 



This fairly general relation can be derived from fac- 
torization properties of partition functions like (10 n]) 
P^ , PH . As a consequence, substituting relation (|13j) with 
7str (x — 2) = — \ in jT2] ) gives the planar bulk dimen- 
sions from the gravity boundary ones: 
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A(°) = U{A) = v(a),V (x) = 1 (4x 2 - 1) 



(14) 



Consider now two sub-stars A, B, issued at the same 
point near a boundary line and avoiding each other. We 
state that their boundary dimensions A on the random 
surface are additive: 



A{AAB) = A{A)+A{B). 



(15) 



Owing to (jl3|), a similar relation, but with a constant 
shift, exists for bulk exponents. On G, two elements in a 
pair {A, B) of walk sets are indeed made non-intersecting 
by gluing a fluctuating patch of random surface between 
them. Each set, Aor B, defines its own independent disc 
as in (j^), with additive b.c.d.'s ([ll 15). So, mutually 
avoiding sets are rendered independent by the fluctuations 
of a random surface. Thus relation ([if]) can be derived 
from factorization properties of exact partition functions 

Consider by contrast two sets A, B of walks which are 
mutually transparent, i.e., A\/ B. In the half-plane R + x 
R, they are independent, and their boundary dimensions 
obey a linear relation 



A< 0) {A\/ B) — A (0) (A) + A (0) (B) 



(16) 



due to the trivial factorization of their partition func- 
tions. On a random surface, their boundary dimensions 
are obtained by inverting ( p^ ) 



A = U-' ( A (0) ) , U 



{x)^-{V24x-+T-l), (17) 



and are not additive. The metric fluctuations indeed cou- 
ple the sets A, B. 

It is clear at this stage that the set of equations above 
is complete. It allows for the calculation of any confor- 
mal dimensions A (S) or A^ ^ (<S) associated with a star 
structure S of the most general type, as in ([!]), involv- 
ing (A,V) operations separated by nested parentheses. 
Any such structure can be systematically reduced: one 
starts from the outermost parenthesis set, and calculates 



b.c.d.'s of operations (A, V) by using (15) for A on the 
random surface G, and (|l6|) for V on the plane R 2 , while 
applying repeatedly the non-linear map U : G — > R 2 
(|l2|), or its inverse U^ 1 ( |l7| ) to transfer to the proper 
space where the boundary dimension is a linear repre- 
sentation of A or V. At the end, one uses ( |l4| ) to recover 
bulk dimensions. 

Brownian-polymer exponents: The single extremity 
scaling dimensions are for a RW or a SAW near a Dirich- 
let boundary in R 2 MM 



>(0) 



(1) = 1, Af (i) 



(18) 



or on G, using (0), A B (1) = U~ l (1) = 1, A P (1) = 



these are the only numerical seeds, i.e., generators, we 
need. 

Stars can include bunches of n copies of transparent 
RW's or m transparent SAW's. Their b.c.d.'s in R 2 
are respectively, by using @ and |ij), A^ } (n) = n 
and Ap^ (m) = |rn, from which the inverse mapping 

( |l7| ) to the random surface yields Ab {n) = U^ 1 {n) 
and Ap (m) = U^ 1 (§?rt) . The star made of L bunches 
i £ {1,...,L} of ni transparent RW's each, and V 
bunches I' £ {1, L'} of me' SAW's, all mutually avoid- 
ing, has planar scaling dimensions owing to ( |l2| , p~4[ ) and 



{m,mi>} = u(a), A<°) { nti mt>} = V ( A 



A{nt,mt>} = 



u- 1 M + 



U- 1 (f m,) 



These exponents are invariant under permutation of the 
bunches of walks. The existence of such a relation has 
been found for RW's in J24|, but with an unspecified U, 
which is here derived from quantum gravity and general- 
ized to SAW's. 

For a copolymer star Sl.l 1 made of L RW's and 
L' SAW's, all mutually avoiding {y£,t,ne = my = 1) , 
A (S l ,l>) =L+\L' gives the c.d. in R 2 



A (0) {S L , L .) = \ 
A<°> {S L , L ,) = i 



(L + |L') (1 + 2L+|L') 



AIL- 
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Because of the star algebra described above 



recovering for L — the SAW exponents |To| ] and for 
L' = the RW non-intersection exponents [|14|. 

Disconnection exponents: Take any walk star A, with 
elementary boundary c.d. x in R 2 . The star S n = {yA) n 
made of n transparent copies of A has b.c.d. nx in R 2 , 
thus U^ 1 {nx) on G. Its bulk c.d. in R 2 is, according to 
d, [l|), x (<S„) = 2V [U- 1 {nx)} , which differs from nx, 
where x = 2V [U^ 1 (x)] is the bulk c.d. in R 2 . Their 
difference is the disconnection exponent, governing the 
conditioned probability Pr ~ j^-x(s n )+nx ^ e un j on 
of n copies does not disconnect the star origin from in- 
finity, within a disc of radius R. 

Multifractal harmonic measure: The harmonic mea- 
sure H {w) of a given set is the probability that a RW 
coming from infinity first hits the set (the absorber) at 
point w. When the set is a RW or a SAW of size R, 
the site average of its moments H n has been shown || 
to be represented by a copolymer star partition func- 
tion of type ([!]) where the absorber avoids a bunch 
of n independent RW's. More precisely J2 W H n (vj) ~ 
Zr {S A n) /Zr (S A 1) , where the absorber S is either 
the two-RW star B V B or the two-SAW star PAP, 
made of two non-intersecting SAW's. We have intro- 
duced the short-hand notation S A n = S A {yB) n de- 
scribing the copolymer star made by the absorber iS 
hit by the bunch {yB) n at the apex only. Owing to 



3 



Eq.(||), we get the scaling ^2 w H n (w) ~ R~ T ( n \ where 
t (n) = (n — 1) D (n) = x (S^n) — x (S A l) defines (an- 
nealed fpl) generalized dimensions D (n) . Our formal- 
ism dl4pq , |l7| ) immediately gives the scaling dimensions 

x (<S A n) = 2V (A (S) + U- 1 (n)) , where A (S) is as 
usual the quantum gravity boundary dimension of the 
absorber S alone. A simple calculation gives r (n) and its 
Legendre transform / (a) + r (n) = an, a — dr (n) /dn 



r(n) =i(n-l) + yi(x/24^+I-5), 



/(«) 



1 

24 



25 
T 



1 



+ 5y--y< (2a -1)- 



(19) 
(20) 



where y = 4A (5) — 1 is the only parameter encoding 
which absorber we consider (which can actually be any 
star tip). For a RW absorber, we have A(B\/B) = 
U- 1 {2) = §, while for a SAW A (P A P) = 2A Pil = 
2U~ 1 (|) = |, thus y = 5 in both cases. The coincidence 
of these two values tells us that in 2D the harmonic mul- 
tifractal spectra / (a) of a random walk or a self-avoiding 
walk are identical. Their Mf spectra associated with walk 
ends j|] however differ, and are obtained using y = 3 for 
a RW end, or y = 2 for a SAW end. 




FIG. 1. Harmonic multifractal dimensions r(n) and spec- 
trum f(a) of a two-dimensional RW or SAW. 

The corresponding universal curves for y = 5 are shown 
in Fig.l: r (n) is half a parabola, and / (a) a hyperbola. 
D (1) = t' (1) = 1 is just Makarov's theorem the di- 
vergence of / at a m in = 5 corresponds to singular needles 
in the absorber, while — r (0) = sup a / (a) = f (3) = | is 
the Hausdorff dimension of the Brownian frontier or of a 



SAW. Thus Mandelbrot's classical conjecture identifying 
the latter two is generalized and proven for the whole 
/ (a) harmonic spectrum. 

Acknowledgements. I thank T.C. Halsey for discus- 
sions and a careful reading of the manuscript. 



[4 
[5 
[6 

\r 

[§ 
[9 
[10 

[11 
[12 

[13 

[14 
[15 

[16 

[17; 

[18 

[19 

[20 

[21 
[22 

[23 
[24 



B.B Mandelbrot, J. Fluid. Mech. 62, 331 (1974). 

H. G.E. Hentschel and I. Procaccia, Physica (Amsterdam) 
8D, 835 (1983). 

U. Frisch and G. Parisi, in Proceedings of the In- 
ternational School of Physics "Enrico Fermi", course 
LXXXVIII, edited by M. Ghil (North-Holland, New 
York, 1985) p. 84. 

T.C. Halsey, M.H. Jensen, L.P. Kadanoff, I. Procaccia, 
and B.I. Shraiman, Phys. Rev. A33, 1141 (1986). 
D. Bernard, K. Gaw§dzki, and A. Kupiainen, Phys. Rev. 
E54, 2564 (1996). 

T.C. Halsey, B. Duplantier, and K. Honda, Phys. Rev. 
Lett. 78, 1719 (1997); J. Stat. Phys. 85, 681 (1996). 
M. Cates and J.M. Deutsch, Phys. Rev. A35, 4907 

(1987) ; B. Duplantier and A. Ludwig, Phys. Rev. Lett. 
66, 247 (1991); C. von Ferber, Nucl. Phys. B490, 511 
(1997). 

M.E. Cates and T.A. Witten, Phys. Rev. Lett. 56, 2497 
(1986); Phys. Rev. A35, 1809 (1987). 
B. Duplantier, Phys. Rev. Lett. 57, 941 (1986); J. Stat. 
Phys. 54, 581 (1989). 

B. Duplantier and H. Saleur, Phys. Rev. Lett. 57, 3179 

(1986) ; H. Saleur, J. Phys. A19 L807 (1986). 

B. Duplantier, Commun. Math. Phys. 117, 279 (1987). 
G.F. Lawler, Intersections of random walks (Boston, 
Birkhauser, 1991). 

B. Duplantier and K.-H. Kwon, Phys. Rev. Lett. 61, 2514 

(1988) . 

B. Duplantier (Saclay preprint SPhT 98/063). 

J.F. Joanny, L. Leibler, and R.C. Ball, J. Chem. Phys. 
81, 4640 (1984); L. Schafer, U. Lehr, and C. Kappeler, 
J. Phys. I 1, 211 (1991). 

C. von Ferber and Y. Holovatch, Europhys. Lett. 39, 31 
(1997); Phys. Rev. E56, 6370 (1997). 

B. Nienhuis, in Phase transitions and critical phenomena, 
Vol. 11, edited by C. Domb and J.L. Lebowitz (Aca- 
demic, London, 1987). 

A. M. Polyakov, Mod. Phys. Lett. A2, 893 (1987); V.G. 
Knizhnik, A.M. Polyakov, and A.B. Zamolodchikov, 
Mod. Phys. Lett. A3, 819 (1988). 

I. K. Rostov and M.L. Mehta, Phys. Lett. B189, 118 

(1987) . 

B. Duplantier and I.K. Rostov, Phys. Rev. Lett. 61, 1436 

(1988) ; Nucl. Phys. B340, 491 (1990). 
B. Duplantier (to be published). 

F. David, Mod. Phys. Lett. A3, 1651 (1988); J. Distler 
and H. Kawai, Nucl. Phys. B321 509 (1988). 

J. Cardy, Nucl. Phys. B240 [FS12], 514 (1984). 

G. F. Lawler and W. Werner (unpublished). 



4 



[25] N.G. Makarov, Proc. London Math. Soc. 51, 369 (1985). 



5 



